x0 = [0.5,10,10,10,10];
options = optimset('MaxFunEvals',1000,'Display','Iter');
[x,fval] = fminsearch(@secretballot,x0,options);
% 
p=(1+tanh(x(1)))/2.0
k=x(2)